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Abstract. The Australian SKA Pathfinder (ASKAP) will be producing 2.2 terabyte 
HI spectral-line cubes for each 8 hours of observation by 2013. Global views of spec- 
tral data cubes are vital for the detection of instrumentation errors, the identification of 
data artefacts and noise characteristics, and the discovery of strange phenomena, unex- 
pected relations, or unknown patterns. We have previously presented the first frame- 
work that can render ASKAP-sized cubes at interactive frame rates. The framework 
provides the user with a real-time interactive volume rendering by combining shared 
and distributed memory architectures, distributed CPUs and graphics processing units 
(GPUs), using the ray-casting algorithm. In this paper we present two main extensions 
of this framework which are: using a multi-panel display system to provide a high res- 
olution rendering output, and the ability to integrate automated data analysis tools into 
the visualization output and to interact with its output in place. 



1. Intoduction 

Upcoming radio observation facilities such as the Australian Square Kilometre Array 
Pathfinder (ASKAP), the MeerKAT Karoo Array Telescope, the Low Frequency Array 
(LOFAR), and ultimately the Square Kilometre Array (SKA) will pose a significant 
challenge for current astronomical data analysis and visualization tools. The expected 
data product sizes (e.g. 2.2 TB per 8 hour of observation for the proposed WALLABY 
all-sky HI survey^) are orders of magnitude larger than astronomers, and existing as- 
tronomy software, are accustomed to dealing with. 



In Hassan et al. (201 1 1, we presented a distributed GPU framework to interactively 
volume render larger-than-memory astronomical data cubes. Throughout this work, we 
demonstrated how volume rendering offers an alternative to standard 2D visualization 
techniques and provided a way to overcome the technological barrier caused by the 
computational requirement of volume rendering for large data cubes. The presented 
framework utilizes a heterogeneous CPU and GPU hardware infrastructure, combining 
shared- and distributed-memory architectures, to yield a scalable volume rendering so- 
lution, capable of volume rendering image cubes larger than a single machine memory 
limit, in real-time and at interactive frame rates. The usage of GPUs as the main pro- 
cessing backbone for the system, and the remote visualization architecture adopted in 



'See 



http://www.atnf.csiro.au/research/WALLABY 
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our design for this framework enables further enhancement for the astronomer's visu- 
alization experience. 

In this paper we present an extension to this framework to provide: better visu- 
ahzation output by integrating external quantitative information with the volume ren- 
dering output (see section [2]) , and high resolution output via multi-panel displays (see 
section [3). 



2. Integrating Duchamp Output 

While the main processing burden is moved to the back-end cluster of GPUs, the ren- 
dering client is free to do some relatively smaller processing to enhance the interactivity 
and the scientific outcomes of the visualization output. The usage of remote visualiza- 
tion enables the viewer application to integrate different graphical primitives with the 



volume rendering output without affecting the rendering processes [see Figure 2 in Has 



san et al. (2011) for details of the framework's main components]. We already utilized 



this ability to present simple primitives such as the colour map, and orientation arrows 



in the original implementation [see Figure 1, 5, and 6 in Hassan et al. (201 1 1]. We have 
extended our work to enable the overlay of a source finder output catalogue on the vol- 
ume rendering output. We use the Duchamp source finder (Whiting 2008), the selected 
source finder for the proposed ASKAP survey pipeline, to demonstrate this integration. 

First, the user loads a Duchamp output file - an ASCII file that contains a set of 
detected astronomical sources represented by their bounding cubes. The viewer appli- 
cation parses the file and extracts the sources' bounding cubes, which are superimposed 
as wireframe boxes on the volume rendering output. The user interaction with the dis- 
played volume rendering output is automatically applied to the displayed source cubes. 
The user can select one of the source cubes to display further available information 
about that source. See Figure [T] for an illustration of this process. The user can also 
show/hide a group of the displayed sources to enable better understanding of the out- 
put. We anticipate other applications of this extension such as comparing multiple 
source finder outputs (by assigning different colour for the source cubes), and overlay- 
ing radio/optical catalogues. 



3. OptlPortal Integeration 

Another benefit of using GPUs is the ability to render high resolution output in a rel- 
atively short time. So, as a further extension of the current viewer implementation, 
we present support for high resolution tiled display systems (such as the OptlPortal 
displa>[^. We used the OptlPortal display facility operated by the Australian Common- 
wealth Scientific and Research Organization (CSIRO) at Marsfield (Sydney, Australia) 
to demonstrate this integeration (see Figure |2]-a). The current CSIRO OptlPortal con- 
sists of 5 X 5 high definition screens (1920 x 1080 pixel) with an overall output size 
of 9600 X 5400 pixels. The viewer uses the back-end GPU cluster to render a static 
image of the current data file with the same orientation information but with a different 
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Figure 1 . Illustration of the process of integerating Duchamp catalogue with the 
volume rendering output. The data cube used is a neutral hydrogen (21cm) observa- 
tion of part of the Ursa Major galaxy cluster, made with the Lovell radio-telescope 
at Jodrell Bank, Manchester. Data courtesy Virginia Kilborn (Swinburne). 



output size (defined by tlie user). Tiie renderingoutput is stored in TIFF format, wliicli 
is converted to tlie pyramidal tiled TIFF formajj 

Figure[2]-b and[2]-c shows an illustration of the importance of high resolution output 
with the HIPASS southern sky cube|^as an example. With the expected increase in the 
output pixel resolution within the different main ASKAP HI surveys (e.g. WALLABY 
- expected output cube size of 6144 x 6144 x 16384), it will be useful to utilize such 
high resolution output facilities to provide the astronomers with a better level of detail 
for both qualitative and quantitative data visualization. 

On the other hand, such facilities may be useful to support collaborative visualiza- 
tion if a suitable data interaction mechanism is provided. We think the current frame- 
work can support an interactive OptlPortal display of the volume rendering output, but 
with a lower frame rate. 



4. Conclusion 

We demonstrate the ability to extend our GPU volume rendering framework to offer 
better visualization outcomes. Two main features were presented: integrating source 
finder output, and utilizing multi-panel displays. Both of these features demonstrate 
two main strong points in the framework design, namely the usage of GPU cluster 
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http://iipimage. sourceforge.net/documentation/iraages/ 



"^The southern sky cube was generated by Russell Jurek (ATNF) from 387 individual cubes. Its dimensions 
are 1721xl721x 1025 with a data size of 12 GB. 
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(b) (c) 

Figure 2. (a) The HIPASS cube in 8000 x 8000 pixel resolution displayed using 
the CSIRO OptlPortal ATNF, Marsfield (9600 x 5400 pixel). Credit: C. J. Fluke, (b) 
A portion of the HIPASS volume rendering with 3000 x 3000 output resolution, (c) 
The same portion of the HIPASS cube with 8000 x 8000 otuput resolution. 



as the main processing back-bone, and the separation between the rendering and the 
result display. We anticipate the ability to further enhance this for integrating more 
quantitative visualization tools and a better utilization of the GPU processing power. 
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